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Abstract 



Strongly coupled positronium, considered in its pseudoscalar sector, is studied 
in the framework of relativistic quantum constraint dynamics. Case's method 
of self-adjoint extension of singular potentials, which avoids explicit introduc- 
tion of regularization cut-offs, is adopted. It is found that, as the coupling 
constant a increases, the bound state spectrum undergoes an abrupt change 
at the critical value a = a c = 1/2. For a > a c , the mass spectrum displays, 
in addition to the existing states for a < a c , a new set of an infinite number 
of bound states concentrated in a narrow band starting at mass W = 0; all 
the states have indefinitely oscillating wave functions near the origin. In the 
limit a — > a c from above, the oscillations disappear and the narrow band of 
low-lying states shrinks to a single massless state with a mass gap with the 
rest of the spectrum. This state has the required properties to represent a 
Goldstone boson and to signal spontaneous breakdown of chiral symmetry. It 
is suggested that the critical coupling constant a c be viewed as a possible can- 
didate for an ultra-violet stable fixed point of QED, with a distinction between 

two phases, joined to each other by a first-order chiral phase transition. 
PACS numbers : 03.65.Pm, ll.10.St, 12.20.-m, 11.30.Rd. 
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I. INTRODUCTION 



The question of a possible existence of an ultra-violet stable fixed point in QED was 
investigated long ago by Gell— Mann and Low [|IJ and developed later by several authors 
0-^]. If such a point were to exist, then the electron mass would be entirely dynamical 
in origin |0-ffl, with a possible spontaneous breakdown of chiral symmetry ||. Although 
perturbation theory calculations do not seem to point to the existence of such a solution, 
quenched lattice QED calculations displayed the existence of a phase transition at the critical 
value a c ~ 0.3 of the coupling constant a, with the occurrence of a spontaneous breakdown 
of chiral symmetry |J. These observations were also confirmed with unquenched lattice 
calculations f7], with a c ~ 0.4, but a vanishing of the Callan-Symanzik function (3 was not 
found there and the question of the validity of QED as a nontrivial consistent theory in the 
continuum limit was raised. 

On the other hand, it is possible, in the continuum theory, to analyze a partial, but 
simpler, aspect of the phase transition problem, namely that of the bound state of strongly 
coupled positronium. If the ground state of the corresponding spectrum, for some value of 
a, were massless, then this would be the signal of a possible spontaneous breakdown of chiral 
symmetry and the prelude of a phase transition in QED. This result would not be, however, 
sufficient by itself to ensure the vanishing of the Callan-Symanzik /3-function, which should 
be shown by independent calculations; only in this case could the consistency of the whole 
procedure be guaranteed. 

We studied, in previous work |||§, the problem of strongly coupled positronium in the 



framework of relativistic quantum constraint dynamics (RQCD) JT0|,|TTJ. This framework 
provides a manifestly covariant three-dimensional description of the internal motion of two- 
body systems and can be shown to be equivalent to a three-dimensional reduction of the 



Bethe— Salpeter equation [T^j. When the approximation of local potentials is made, the 
corresponding wave equations can be analyzed rather easily and in many cases analytic 
solutions can be obtained. It was found that the Todorov form of the electromagnetic two- 
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body potential, first introduced in the quasipotential approach leads to the existence 
of a critical value of the coupling constant, with a c = 1/2. For a > a c , the potential 
becomes too singular and needs some regularization, provided by a cut-off radius tq in the 
Coulomb interaction. In the regularized theory the bound state mass spectrum displays, 
for a ^ 1/2, a rapid fall of the ground state mass to values close to zero, thus indicating a 
drastic modification of the qualitative features of the bound system and presumably of the 
theory itself. The above procedure requires, however, a numerical treatment of the equations 
and makes it difficult to ascertain the existence of solutions for vanishing values of tq. 

The presence of a cut-off radius in the Coulomb interaction needs, in general, an adequate 
interpretation. As long as one studies the behavior of charged particles in supercritical 
Coulomb fields of heavy nuclei, the finite size of the latter naturally regularizes the Coulomb 
interaction at short distances. If, on the other hand, it is the strong interaction of pointlike 
particles that is considered, as in the case of positronium in QED with strong Coulomb 
coupling, then the meaning of r remains unclear. 

A similar problem also occurs with the Bethe-Salpeter equation in the ladder approxi- 



mation. It has been shown [H] that, when the coupling constant a is larger than a critical 
value a c (~ 7r/3 in the Landau gauge and ~ tt/4 in the Feynman gauge), the theory un- 
dergoes spontaneous breakdown of chiral symmetry. However, for a ^ a c , the treatment 
and resolution of the equation necessitate the use of an ultra-violet cut-off A. While the 
introduction of the latter can naturally be justified in QCD as being an approximate way 



of parametrizing the asymptotic freedom of the theory fl5f , it has not received a simple 



interpretation in QED. In this respect, Miransky et al. (T^j suggested that, for a > a c , QED 



undergoes an additional charge renormalization that absorbs the infinities of the pointlike 
limit; the renormalized charge remains equal to a c and then might be identified with the 



ultra-violet fixed point of QED. It was also pointed out in this connection [p~7| , [T8[1 that, 
for a = a c , because of the new renormalization of a, the fermion composite operator ipij) 
acquires the dimension 2 instead of 3, and thus allows for the presence of renormalizable 
four-fermion interactions. 



In quantum mechanics, there exists an alternative method for dealing with singular 
interactions of pointlike particles, without the need of introducing cut-offs: this is the self- 
adjoint extension method, discussed a long time ago for singular potentials of the type l/r n 
(n > 2) by Case |Tj|. In his classic paper, Case showed that all self-adjoint extensions of 



the Klein-Gordon-Coulomb (or Dirac-Coulomb) problem can be parametrized by a single 
constant B when the interaction becomes singular. 

While Case's method is only of academic interest for problems concerning supercritical 
Coulomb fields of heavy nuclei, the size of the latter providing a natural short distance cut- 
off, it reveals its full power in the present problem of strong Coulomb interaction of pointlike 
particles. The self-adjoint extension parameter B can be interpreted as parametrizing the 
short distance behavior of the interaction and its choice amounts to fixing the energy of one 
of the bound states (the values of the masses of the constituent particles of the bound states 
and of the coupling constant a being already fixed) and calculating the other bound state 
energies with respect to this one, without making explicit cut-offs appear. 

In the case of ^-function interactions in two and three space dimensions, it was shown 
PD|| that the self-adjoint extension method provides the renormalized version of the theory, 
when the cut-off of the regularized theory is removed and a corresponding renormalization 
of the coupling constant is performed. It is then natural to expect from the same method of 
approach, applied now to the 1/r 2 singularity, to also provide the finite renormalized version 
of the theory, provided one of the bound state energies is fixed. 

Motivated by these results, we have investigated with Case's method of self-adjoint exten- 
sion the problem of strong Coulomb coupling in positronium-like systems. The relativistic 
wave equations of constraint theory lead for the relative motion in 1 5'o states to a final three- 
dimensional equation which is very similar in form to the Klein-Gordon (KG) equation ||, 
and therefore Case's method can be readily applied to it. 

Our main results are the following. We find that the system undergoes a first-order chiral 
phase transition at the critical value a = a c = 1/2. While the ground state mass for a < 1/2 
can be continued to the domain a > 1/2 and remains different from zero, a new set of an 



infinite number of states, concentrated in a finite domain of mass with accumulation at the 
value zero, appears, the zero mass state representing the new ground state of the system. 
This result occurs for any fixed value of Case's constant B. All the states have indefinitely 
oscillating wave functions near the origin. While tachyonic solutions formally exist, they are 
ruled out from the spectrum by the self-adjointness condition, and therefore the zero mass 
state remains the physical ground state of the spectrum. 

In the limit a — > a c from above, the short-distance oscillations disappear from the wave 
functions, and the states accumulated around the zero mass solution shrink to a single 
massless state with a definite mass gap with the rest of the spectrum. It turns out that the 
latter state has the required properties to represent a Goldstone boson and hence to signal 
a spontaneous breakdown of chiral symmetry. The fact that for a > a c a sensible theory, 
with finite and nonvanishing couplings to the observable currents, can be defined only for 
a = a c + strongly suggests the identification of a c with an ultra-violet stable fixed point of 
QED, with the distinction between two phases, governed by a c — and a c + 0, respectively, 
and joined to each other by a first-order chiral phase transition. 

The paper is organized as follows. Section II discusses Case's method for Klein-Gordon 
particles in an external Coulomb field and is included to make the paper self-contained. 
Section III is devoted to the study, in the framework of RQCD, of the strongly coupled 
positronium spectrum in its pseudoscalar sector. Section IV deals with the question of an 
eventual appearance of tachyonic states in the spectrum of states. In Sec. V, the limit 
a — > a c + is considered and the presence of a Goldstone boson established. Summary 
and discussion of results follow in Sec. VI. In the appendix, some results of RQCD are 
summarized. 



II. CASE'S METHOD FOR THE KLEIN-GORDON EQUATION 

For the sake of completeness and for easy comparison with RQCD, we discuss here Case's 
method |TP[ for a Klein-Gordon particle of charge e, mass M and energy E in a Coulomb 
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potential V(r) = —aZ/r. We have a = 1/137 and fi = c = 1 in our system of units. The 
radial KG equation for s-states is : 

2 E a Z a 2 Z 2 



u" + 



i? — M + 



+ 



u = 



(2.1) 



From Eq. (2.1) one derives the orthogonality condition for two solutions u\ and u 2 of energy 
E 1 and £ 2 (£i ^ E 2 ) : 

2aZ s 



(mi,m 2 ^ 



KG 



- (u 2 'ui - UiU 2 ) 



U\U2 dr 



. 



(2.2) 



Note further that the KG norm is given by 



K u , u )kg 



aZ\ , 



E H u l dr . 



(2.3) 



For aZ > 1/2, the square integrable solutions of Eq. (2.1) (vanishing at r = oo) are 
given by 



u{E,r) = cW k ^{p) , 



(2.4) 



where c is a normalization constant, Wk )fl is the Whittaker function || and where 



p = 2(M 2 -E 2 )ir 



(2.5) 



EaZ 



(M 2 -E 2 )-2 



(2.6) 



ji = i A , A = (a Z 



(2.7) 



Using the formula pi 



r(|-/i-A:) r(| + ^-/c 



(2.8) 



one finds : 
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limw ~ p2 cos (8 + A hip) , (2.9) 

with 

r(-2iA) 

It is easy to verify that the general behavior of the solution of Eq. (2.1) close to r = is 
given by |19| : 

limu ~ cos(Aln(Mr) + B) , (2.11) 

where B is an arbitrary constant. Comparison of Eqs. (2.9) and (2.11) yields the relation : 

argT(l-2zA) + | - argT (- - i A - A;) + Aln2 

(M 2 - E 2 ) 3 

+ Aln^ -^-L- = B + nir, (2.12) 



where n is an arbitrary integer. Furthermore, it can be shown |19| that keeping the same 
value of B for all states guarantees the orthogonality condition (2.2). Therefore, choosing the 
value of B provides a self-adjoint extension of the KG equation and allows us to obtain the 
corresponding spectrum by solving Eq. (2.12) for E. Notice that this equation is invariant 
under adding to B any multiple of ir and, therefore, it is sufficient to consider the values of 
B in the interval [0, n]. 

It is of interest to investigate the small A behaviour of (2.12). Using the formula |E| : 

argT( x + iy ) — y ib (x) + f — - — — arctan — - — ] , (2-13) 

n = o \ x + n x + n J 

where 

T' (r ) 

one finds from Eq. (2.12) : 



-2AV(1) + \ + AV(~-A0 + Aln2y / 1-^ + 0(\ 2 ) = B . (2.15) 
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(Without loss of generality, we can drop here the nir term for small values of A.) For this 

71 

equation to be satisfied (for B ^ — ), as A — > 0, it is necessary that 

fimA^-AO = £-| , (2.16) 

or that 

limk = p+^, (2.17) 

where p is any non negative integer. Using Eq. (2.6), it is easy to see that one recovers the 
usual (£ = 0) spectrum of the KG equation for aZ tending to 1/2 from below. This means 

1 7T 

that the spectrum is continuous through aZ = — for any value of B ^ — . However, the 
slope of the energy curves (in A or aZ) are not continuous. For small values of A, one can 
show that the eigenvalues correspond to the following behavior : 

7T 

On the other hand, for B = — , another state with no correspondence with the spectrum 
for aZ < - appears. Its energy is given, from Eq. (2.15), by the equation 



1 E 2 
iK- - A;) +ln2y 1 — — - 2^(1) = , (2.19) 

and has the value E/M ~ -0.049 . 

We now turn to the question of a possible instability of the system, described by the 
KG equation (2.1) with a pointlike attractive Coulomb interaction, with respect to the 
spontaneous pair creation. This would correspond to the existence of a critical value of 
Z = Z c such that E (Z c ) = — M, or, according to Eq. (2.6), k = — oo. To see whether such a 
solution of Eq. (2.12) exists for small values of A, it is advantageous to transform the third 



term of Eq. (2.12) by using the following formula [23 



Im In T (x + iy) = ar g T (x + iy) = (^x — — ^ arctan — + y | In x 



+ -ln 

2 



l+(- 

x 



y x 2 



with 

x y y x 

a = ~2~, 2 ai 2 6l ' b = 2 01 H 2 61 ' 

ar + y x z + y z x z + y x z + y z 

ai = x2 - y \ , bl = — zxy__ (2 21) 

(x 2 + y 2 ) {x 2 + y 2 ) 

One finds, using Eq. (2.12) for large values of \k\ and Euler's constant 7 = — = 
0.57721..., that the critical value of Z c is given by 

2A c7 + A c ln(2aZ c ) + O (A 2 ) = B - |, (2.22) 

where A c = (a 2 Z^ — 1/4^ ^ . Note the cancellation between the two terms containing 
In (M 2 - E 2 ) in Eq. (2.12). From Eq. (2.22), we see that a critical value of aZ = aZ c > 1/2 

7T 7T 

can only occur for B 7^ — . Also note that A c is exactly if B — — , and will be larger and 
larger if B is increasing. There is no solution of Eq. (2.22) (at least for small values of A c ) 

7T 

for B < — , as the left-hand side should be positive. The above comments are illustrated in 
Figs. 1 and 2 which show respectively how the energy of the lower bound states vary with A 

71 71 

for B = 1.56 < — and B = 1.58 > — respectively. One can see that no instability will occur 

71 71 

(i.e., no state with E = —M) for A < 2 in the first case (B < — ), while, for B > — , there 
is a state with energy E = —M for A ~ 0. Furthermore, the solution given by Eq. (2.19) 
corresponds (for small A) to the boundary case B = n/2 between the two domains B < n/2 
and B > n/2. 

For the KG equation, the sign of the KG norm unambiguously distinguishes between 
particle and possible antiparticle bound states. We have calculated numerically the norm of 
the states of E ~ — M wave functions and established that they are indeed particle states. 
We further show in Fig. 3 the radial wave function corresponding to the lowest state for 
aZ = 1 and B = 1.58 (the dot-dashed curve in Fig. 2). It, of course, displays an infinite 
number of nodes consistent with the behaviour (2.9), but the outermost of these nodes 
occurs at a distance representing a rather small fraction (here of the order of one tenth) of 
the "radius" of the state. Globally, the wave function has a "nodeless" structure, typical of 
a ground state wavefunction. 
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We thus find that the stability property with respect to pair creation of a Coulomb 
source with vanishing radius and charge Z > 137/2 crucially depends on the choice of the 
self-adjoint extension characterized by the constant B. In principle, this constant should be 
determined by comparison with experimental data. However, the case of physical interest 
corresponds here to a finite radius R of the source and has been extensively discussed in 
the literature ^4|. We shall not pursue further the investigation of this academic problem. 
Let us simply note that the Coulomb field of a charge with radius R satisfying MR « 1 
is supercritical for otZ > aZ c = 1/2 f8j, so that the choice B > n/2 would be mandatory to 
reproduce finite radii results. 



III. STRONGLY COUPLED POSITRONIUM 



We now consider the problem of two particles of equal mass m and opposite charges with 
spin one-half, in mutual electromagnetic interaction. Within RQCD the 1 Sq states of the 



system with the Todorov choice of the interaction |T3[ are described in the cm. frame by 
the radial equation ||[| (see also the appendix) : 

d 2 



dr 



2 - TTl w 



(e w - A(r)Y 







(3.1) 



with 



A(r) 



a: 
r 



m w 



W 2 - 2m 2 
2W 

2 

m 

W ' 



(3.2) 



(3.3) 



(3.4) 



W being the cm. energy of the two-body system. Equation (3.1) correctly describes the 
physical positronium 1 S'o energy levels to order a 4 [|IO| , |rT|| , with a being the fine structure 
constant. Our interest is in the solution of Eq. (3.1) for large arbitrary values of a (strongly 
coupled positronium). From now on, we shall consider a as a free parameter. 
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As Eq. (3.1) is very similar to the Klein-Gordon equation (2.1), Case's method can be 
readily applied to it to obtain its solutions in a similar way. From Eq. (3.1) we find, for two 
square integrable solutions <fi and <p 2 of energy W\ and W 2 (W\ ^ W 2 ), the orthogonality 
condition : 



(^1,^2) = 2 - <Pl<P2) 



= 0, 



where the scalar product is now defined by : 

'W 1 + W 2 



(<Pl,<P2) = / 
JO 



A(r) 1 + 



2m 2 



4 V W 1 W 2i 

The conserved norm iV corresponding to Eq. (3.6) is given by : 



y?iy?2 dr . 



N 



= r 

Jo 



W . . . / 2m 2 



ip 2 dr . 



For a > 1/2, the solutions of Eq. (3.1), vanishing at r = 00, are given by : 



with 



k = 



OLEw 



a{W 2 - 2m 2 ) 



{m 2 w -e 2 w f /2 W(4mi-WY /2 ' 



i A , A 



1\ 1/2 



a "4 



(3.5) 



(3.6) 



(3.7) 



(3.8) 



(3.9) 



(3.10) 



p = 2iTr , 



(3.11) 



K = (m 2 w 



eV)' 2 = \{Am 2 -W 2 ) 



1/2 



(3.12) 



As r — > 0, the wave function </?(r) exhibits the behavior described in Eqs. (2.8), (2.9) 
and (2.11) (with k replaced by k). Choosing the same value of B for all the wave functions 
guarantees the validity of the orthogonality condition (3.5) for W\ 7^ W 2 - The corresponding 
energy spectrum is then given by a relation similar to Eq. (2.12) : 
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ar g r(l-2iA) + | - aigT^-iX-k^j + Aln^— ) Z? ;/t . (3.13) 



m 



It is easy to see, using Eqs. (2.13) and (2.14), that, for small values of A, there are 
solutions of Eq. (3.13) corresponding to 

lim k = p+ - , (3.14) 

where p is any non negative integer, meaning that the a < - spectrum || is extended 
continuously through a — -. 

On the other hand, Eq. (3.13) possesses a new set of solutions with energies located 
in the positive vicinity of zero. To exhibit them, we consider small values of A and large 
negative values of k, corresponding to small and positive values of W (see Eqs. (3.9) and 
(3.12)). Using Eqs. (2.20) and (2.21), we can rewrite Eq. (3.13) in this case as 

-2A + Xln(--k) + Alnf— ) + O ( X 2 ) = B - - + nvr , (3.15) 
2 \ m J v ' 2 

or, using Eqs. (3.9) and (3.12), as 

-2A - Alnf— o\ + A In 2 + O (\ 2 ) = B - - + nir . (3.16) 

\m J v ' 2 

There are thus solutions of vanishing W as A — » : 

-Aln — a) -> B + mr , (3.17) 

\m J 2 



W „ a -x p -(B+mr-7r/2)/\ 
m 



where n = 0, 1, 2,... for B > tt/2 and n = 1, 2, 3,... for I? < ir/2 (the reason being that the 
right-hand side of Eq. (3.17) should be positive for small values of W). 

The absence of these new types of solution in the Klein-Gordon equation case is due 
essentially to the fact that in Eq. (2.12) the term ln((M 2 — E 2 ) 1 ' 2 /M) cancels, in the limit 
E — > —M, the contribution coming from the large negative values of k; no such cancellation 
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exists in the present case, for the equivalent parameter K [Eq. (3.12)] remains finite and 
different from zero in the limit W — > 0. 

Equation (3.18) shows that, for a fixed value of B and sufficiently small A, there is an 
infinite number of states concentrated in a narrow band of energy and having the zero energy 
as an accumulation point. Clearly, the zero energy state is the new ground state of the bound 
state spectrum. The energy spectrum given by Eq. (3.13) is plotted in Figs. 4 and 5 and 
illustrates the above comments again for values of B < ir/2 (Fig. 4) and B > ir/2 (Fig. 5). 
Only a few states corresponding to Eq. (3.18) are shown. 



IV. THE FATE OF TACHYONS 

To complete the analysis of the energy spectrum, we discuss in this section the question 
of the status of tachyonic states. The presence of such states in the spectrum of physical 
states would invalidate the conclusions of Sec. Ill, since the zero energy state would no 
longer represent the ground state of the spectrum, the latter then displaying rather the 
characteristics of an unstable vacuum. 

Tachyons are characterized by negative values of W 2 , i.e., by imaginary values of W. 
Inspection of Eqs. (3.9) and (3.12) shows that in this case k becomes imaginary, while K 
remains real. 

Equation (3.1) still has the Whittaker functions as normalizable solutions. We 

define 

k = if] , (4.1) 

and examine the behavior of the wave function near the origin; it is given by Eq. (2.8) : 

r(-2iA) r(2iA) i / 

^ * tvi r, — p 1+ + tvi rr, p"~ iX ■ (4-2) 



In order to display the main qualitative features of this type of solution, we consider, in 
the following, large values of \q\, while keeping A small. Then, the functions T(- ± iX — irj) 
can be expanded in terms of iX : 
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r(i±zA-z»j) ~ r(~-ir ] )e ± ( iXln \ r >\ + <v)^/2) ; (4.3) 
[e(r?) is the sign of rj,] leading in turn to the following behavior of W~ : 

w~ k J P ) - f( r 1 -y { r(-2<A) e iX ln N + e ^/ 2 ^ 

+ T(2zA) e~ iX ln M - e M A7r / 2 pH* ) . (4.4) 



The wave function if has the same formal behavior as in Eqs. (2.9) or (2.11), with /3 
defined as 

f3 = argT(-2iA) + Aln|r/| -ie{r])X^ . (4.5) 

Notice that (3, and hence B, is complex. 

In order to study the orthogonality conditions for these states, we observe that when 
tachyons exist, there is a doubling of states : for each "eigenvalue" irj, the value — irj is also 
a solution. The rules of constructing scalar products for such states have been studied in 
the two papers of Ref. P5] . For a state with complex "eigenvalue" irj one also considers 



the state with the complex conjugate "eigenvalue" —irj, called the associated vector. While 
the norm of a state with complex eigenvalue is zero, its scalar product with its associated 



vector is nonzero f25| . Therefore, the following rule should be adopted for the choice of 
admissible states. Define admissible states as those corresponding to a definite sign of rj in 
the "eigenvalues" irj. Then, the corresponding adjoint states in the scalar product are the 
associated vectors (with "eigenvalues" —ir/); these appear there complex conjugated. 

The analysis of Case can then be repeated for the tachyonic states. Once the same 
complex value of B is chosen for the tachyonic solutions, then these will satisfy among 
themselves the orthogonality conditions (3.5) and might constitute admissible states. [Notice 
that associated vectors cannot be considered as admissible states, for, according to Eq. (4.5), 
they would not have the same value of B; they appear only as the adjoints of admissible 
states.] However, the orthogonality condition (3.5) fails when considered between a tachyonic 
state and a "normal" state (with W 2 > and W > 0), because the coefficient B cannot be 
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chosen the same for both types of solution (it must be complex for the former and real for 
the latter). This is the reason why tachyons are rejected from the Hilbert space of physical 
states. 

Therefore, we end up with the conclusion that the physical spectrum is free of tachyons 
and its ground state is the zero energy state. 

V. THE LIMIT a ->■ - + e AND SPONTANEOUS CHIRAL SYMMETRY 

BREAKING 

The occurrence of a zero mass ground state, with the quantum numbers of a pseudoscalar 
boson, in the bound state spectrum is suggestive of spontaneous chiral symmetry breaking. 
However, the oscillatory nature of the wave functions near the origin makes it difficult 
to define their couplings to the axial vector current. Furthermore, the accumulation of 
an infinite number of states around the zero mass state [Eq. (3.18)] does not allow the 
disentanglement of the ground state from the rest of the states of its neighborhood. 

Within these circumstances, the limiting value a = - + e (e = +0) plays a particular role 
for a physical interpretation of the theory. In this limit, the oscillating behaviors of the type 
(2.9) or (2.11) disappear from the wave functions and, according to Eq. (3.18), all the low 
mass states shrink to a single state with zero mass. A definite mass gap appears between 
the zero mass ground state and the other massive states of the spectrum. 

To study in more detail the properties of the corresponding wave functions, it is preferable 
to reanalyze Eq. (3.1) for the particular value a — - + e. The normalizable solutions are 
the Whittaker functions W~ , which behave near the origin as 

W~ kfi (p) ~ (ln2p + ^(-fc + l)-2lKl)) , (5-1) 

the various parameters and variables being already defined in Eqs. (3.9)-(3.12), while the 
general behavior of the solutions of Eq. (3.1) is : 

ip ~ r 1/2 (aln(mr) + b) , (5.2) 
16 



where a and b are constants. 

The orthogonality condition (3.5) requires from the admissible normalizable solutions to 
satisfy the condition 

b 



a 



A , (5.3) 



A being the same (arbitrary) constant for all the solutions. Eq. (5.1) then yields the 
eigenvalue equation : 

2K ~ 1 

ln( ) + ^(-k + -) - 2^(1) = A. (5.4) 

m 2 

In order that the corresponding solutions be the limits of those found for a > 1/2 [Eqs. 
(3.14) and (3.18)], it is necessary that the constant A equal +oo. One therefore finds again 
the solutions (3.14), as well as the additional single solution W = 0, corresponding to 
k — > — oo. 

The infinite value of A means that for the solutions of the type (3.14) the logarithmic piece 
in Eq. (5.2) is absent and the corresponding wave functions behave as r 1 ' 2 near the origin. 
This is not true for the solution corresponding to W = 0, for the limit k — > — oo cannot 
be straightforwardly taken in the Whittaker function W^ Q . To analyze more accurately the 
properties of the corresponding wave function, we consider small values of W (W = W — 0), 
make in Eq. (3.1) the change of variable 

m2r k\ 
V = Wo' (5 - 5) 

and keep only leading terms, it being understood that the limit W — > should be taken at 

the end of calculations of physiacal quantities. Equation (3.1) becomes at leading order in 

Wo : 

^ ~ \* + 4?^ = ° ' (5 - 6) 



the solution of which is : 



1 /2 \ 

tpo = corny* e-^ *(-, 1, Ay 1 ' 2 ) , (5.7) 
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where cq is a dimensionless normalization constant and \I/ (also denoted by U in the liter- 
ature) is the confluent hypergeometric function which behaves at infinity with a power law 
pl| ; its dominant behavior near the origin is given by 

lim*(~ l,x) = ( lux + 2^(1)) . (5.8) 

The solution (5.7) is a nodeless function and represents the ground state of the spectrum. 
[The function $(a,c;a;) does not have positive zeros for a and c real and either a > or 
a — c+ 1 > 0. Also, it can be checked that the formal tachyonic solutions, found in Sec. IV, 
are absent in the present case.] It has a distribution-like behavior, due to the fact that it is 
defined by the limiting procedure Wq — > : taking the limit Wq — *> 0, while keeping r fixed, 
shows that the wave function is actually peaked at the origin (recall that the complete wave 
function is <f/r). 

To compute the normalization constant of this wave function (as well as of the others), 
it is necessary to reconstruct the whole sixteen-component spinor wave function and to use 
its relationship with the Bethe— Salpeter wave function, which ultimately fixes the normal- 
ization coefficients. Some details of these calculations can be found in the appendix . One 
finds for the norm N of Eq. (3.17) the expression : 

W 2 

N = — , (5.9) 

where W is the mass of the bound state. 

As to the massless state, after using the change of variable (5.5) in Eq. (3.7), we find 
that the dominant contribution for small W comes from the third term in the integral. This 
leads to the following behavior of the normalization constant Co of the wave function po(y) 
[Eq. (5.7)] for small W : 

c ~ Wljm 2 . (5.10) 

We now turn to the calculation of the coupling constants of the bound states to the axial 
vector current. These are defined as 



<0|j*(0)|P> = P/j,F , (5.11) 

where P is the four-momentum of the pseudoscalar state and is the renormalized axial 
vector current. In general, in the absence of anomalies, the axial vector current undergoes 



only a finite multiplicative renormalization by radiative corrections p6[ . This feature is, 
however, the result of compensating contributions from propagator and vertex renormaliza- 
tions (with factors Z2 and Z^ 1 , respectively). More explicitly, one has 



3n5 - Z 2Z A l j^ 5 , 



(5.12) 



where is the unrenormalized current, with Z^Z A finite. 

From operator product expansion J2TJ and renormalization group analysis ]28j one finds: 



< 0\jg\P > = -Z A T no75 <p BS (x) 



(5.13) 



where <Pbs is the Bethe— Salpeter wave function. Using the relationship of the constraint 
theory wave function ip with the Bethe— Salpeter wave function <Pbs [Eqs. (A6)-(A8)] one 
obtains (in the cm. frame) : 

2a\~ 1/2 



<0\]£\P> =-Z A T no75 1 + 



Wr 



(5.14) 



2a \ - 1/2 

\\\r = 2Z A (l + —j Trip 4 (r) 



(5.15) 



The renormalization constant Z A should render finite the physical coupling constants P. 
Let A be the ultra-violet cut-off (in momentum space) of the four-dimensional theory and let 
r be the corresponding short-distance cut-off (in rr-space) of the three-dimensional theory. 
We shall admit the weak relation 



lim To (A) 

A— >oo 







(5.16) 



and shall transpose several known qualitative results of the four-dimensional theory into the 
three-dimensional one. When tq 7^ 0, the electron has a bare mass mo(ro) that vanishes 
with r with some power v : 
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lim m (r ) ~ mOmroY , v > . (5-17) 

For m 7^ 0, the Ward identities of the axial vector current imply that the Goldstone boson 
acquires a mass Wo, which behaves in terms of mo as [|5|. |T5*| . [29| 

lim W^(r ) ~ m (r )m . (5.18) 

We shall assume that v < 2; because of Eq. (5.18), it is only in this case (including eventually 
the limiting case v = 2) that the wave function (5.7) can be consistently defined. As a matter 
of rough comparison, the analog of v in the four-dimensional theory, calculated at the two 
loop level, is equal to 3a /(2tt) + (3/4)(a/(2vr)) 2 ~ 0.25 for a = 1/2 g|. 

We designate by W\ and F\ the mass and coupling constant of a massive state of the 
bound state spectrum and by W and F the similar quantities of the ground state (the 
Goldstone boson) when r$ ^ 0. Taking into account the behaviors of the corresponding 
wave functions near the origin [(fix ~ mci(mr) 1 ^ 2 , (p ~ mco(m 2 r /Wo) 1 / 2 hi(m 2 ro/Wo)}, Eq. 
(5.15) leads for the two cases to the following equations : 

8m 3 Co\n(m 2 ro/Wo) 
F ° = ZA W 2 iWoroY/ 2 ■ (5 ' 20) 

In Eq. (5.19), C\ is a nonvanishing normalization constant in the limit r — > 0, while the 
bahavior of Co in Eq. (5.20) is given by Eq. (5.10). In order to maintain finite the value of 
F in the limit r —>■ 0, we must have : 

Z A ~ (lyoro) 172 / 111 ^^)" 1 ~ (mro)^ / \n{mro)~ l . (5.21) 

Replacing in Eq. (5.19) yields F\ = 0. Therefore we obtain, when r — > 0, the following 
behaviors of the coupling constants : 

F ^0, F 1 = 0. (5.22) 
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These are precisely the complementary conditions for having spontaneous breakdown of 
chiral symmetry : only the Goldstone boson couples to the axial vector current. 

The behavior of Za when r$ —>■ 0, given by Eq. (5.21), is in qualitative agreement with 
its behavior in the four- dimensional theory : in the Feynman gauge, Z2, and hence Za, 
vanishes when the ultra-violet cut-off A goes to infinity || . 

We can also calculate the matrix elements of the divergence operator <9 M j M s; actually, this 
should only lead to a check of the covariance property of the formalism. We find : 

<0\d»j*\P> =-iW 2 F 



-Z A 2im(l + — ^-) 1 Trj 5 i/j(r) 

-Z A ^m(l + ^)- 1 Trfa(r) 

Wr r _>o 

— Zjfiim— 

r r ^o 



(5.23) 



which yields back Eq. (5.15). [In obtaining the above results, the wave equations of ip have 
been used; similar calculations can be found in the second paper of Ref . [|II| .] According to 
Eqs. (5.22), for all states of the bound state spectrum, the matrix elements (5.23) vanish, 
thus ensuring axial vector current conservation. 

We end up with the conclusion that, in the limit a = - + e, the theory displays the 
features of spontaneous chiral symmetry breaking : i) the presence of a zero mass ground 
state in the spectrum, with a mass gap with the rest of the bound states; ii) a nonvanishing 
coupling of the Goldstone boson to the axial vector current with a decoupling of the massive 
states from the latter. 

For completeness, let us also describe the situation that results from the limit a = l^~ t - 
For a < — , the positronium spectrum has the usual structure ||, the same as for small a, 
without massless bound states, and the limit a — > - from below does not introduce any 
qualitative changes, the ground state remaining massive. In this case the renormalization 
constant Za is determined from the finiteness of Fi [Eq. (5.19)]. The latter should be 
different from zero; otherwise, the axial vector current would be conserved and, in the 
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absence of a Goldstone boson, chiral symmetry would be realized through its normal mode, 
implying a parity doubling of degenerate states; this is not realized in the bound state 
spectrum; furthermore, in QED, with only electrons as massive fermions, this also is not 
possible. We therefore conclude that for a = - — e, the axial vector current is not conserved 
and chiral symmetry is explicitly broken by the electron mass. 



VI. SUMMARY AND DISCUSSION 

We applied Case's method of self-adjoint extension of singular potentials, to the study 
of strongly coupled positronium in its pseudoscalar sector in the framework of relativistic 
quantum constraint dynamics. We found that, as the coupling constant a increases, the 
bound state spectrum undergoes, at the critical value a = a c = 1/2, an abrupt qualitative 
change. For a > a c , the mass spectrum displays, in addition to the existing states for 
a < a c , a new set of an infinite number of bound states, concentrated in a narrow band 
starting at mass W = 0. The bound states have indefinitely oscillating wave functions near 
the origin. 

In the limit a — > a c from above, the short- distance oscillations disappear and the states 
accumulated around the zero mass state, shrink to a single massless state, representing the 
ground state of the spectrum, with a definite mass gap with the rest of the states. This state 
has the required properties to represent a Goldstone boson and hence it signals a transition 
to a new phase where chiral symmetry is spontaneously broken. This is also an expected 
possibility from the existence of an ultra-violet stable fixed point in QED and is therefore 
suggestive of an identification of the critical value a c with the Gell— Mann— Low eigenvalue 
ceo- The fact that the two boundary values a c — e and a c + e correspond to different phases, 
the former governing a phase where chiral symmetry is broken by the electron mass term and 
the latter governing a phase where chiral symmetry is spontaneously broken, necessitates 
the introduction of a similar distinction for ao, with boundary values ao — e and ao + e, with 
the search for the relevant domains in the theory. 
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In the model potential we were considering, the contribution of the one-photon exchange 
diagram, besides kinematic factors, is represented in the three-dimensional theory by a/r, 
which means that no distinction was made between the large- and short- distance behaviors 
of the effective charge (the physical coupling constant a being, in general, determined from 
the large distance behavior of the photon propagator). In order to be able to determine 
the domains of each of the above phases with respect to the values of the physical cou- 
pling constant, we define the effective charge in four-dimensional momentum space |!H3| as 
a e ff(—q 2 /m 2 ) = ad(—q 2 /m 2 ,a); d is the Lorentz invariant part of the transverse part of 
the photon propagator multiplied by q 2 ; a is the physical coupling constant, measured at 
large distances : a e //(0) = a. The asymptotic value of a e ff is a [a e //(oo) = a ]> at which 
value the Gell— Mann— Low function ip vanishes. Because of the positivity of the photon 
two-point spectral function, one has in general the inequality a e ff(—q 2 /m 2 ) < a for q 2 < 
0, which implies in particular that a < «o- 

Another particular value of a is provided by the zero of the Callan-Symanzik function 
(3, which we denote by a 1; with «i < a , satisfying aid as (l, a\) = a , where d as is the 
asymptotic part of d [[J. If radiative corrections are estimated to be of the order of a /it in 
general, then, for a ~ a = 1/2, they are of the order of 20% and a\ should be of the order 
of 0.4. It was shown by Adler |§ that, according to the ways of summing diagrams, either 
cti or a are essential singularities for the corresponding defining functions ((3 or ip). The 
value «i appears then as a natural separation point between two subdomains in the domain 
of variation < a < a^. The following scheme might provide a possible description of the 
conditions of occurrences of each of the phases mentioned above. 

When < a < ai, the asymptotic behavior of the photon propagator is governed by 
a — e and we are in the phase where chiral symmetry is broken by the electron mass. In 
this phase, the renormalized vertex function mT (defined in Ref. 0), corresponding to the 
divergence of the axial vector current, is different from zero and the axial vector current is 
not conserved (in the absence of anomalies). The structure of the positronium spectrum is 
the same as for small values of a, with a massive ground state. 
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When a jumps from a\ — e to a\ + e, this induces, through the relationship between the 
/3 and i[) functions, a similar jump of the bare coupling constant from — e to a® + e and 

— 5 

we enter in the phase where chiral symmetry is spontaneously broken. In this phase mT 
is identically zero and the axial vector current is conserved. The positronium spectrum has 
now, in addition to the existing states for a < a\, a massless ground state. 

When a > «o, because of the positivity condition already mentioned, the unitarity of 
the theory breaks dowm. 



Our conclusions also join those obtained by Miransky et al. |L6[ from the Bethe— Salpeter 
equation, who conjectured that the critical value a c should be identified with the 
Gell— Mann— Low eigenvalue. The difference in the numerical values of a c found in the two 
approaches (a c = 1/2 here and a c = 7r/4 in the ladder approximation of the Bethe— Salpeter 
equation in the Feynman gauge) is presumably related to the different approximations used 
in the kernels of the bound state wave equations. The Todorov potential, used in the present 
approach, takes into account multi-photon exchange diagrams and correctly reproduces the 



physical positronium and muonium spectra to order a 4 |10|JTT[; this is not the case for the 



ladder approximation of the Bethe— Salpeter equation in covariant gauges. 

In the course of the present analysis, the effects of anomalies were ignored; these are 



known to modify the Ward identities of the axial vector current f3C|. However, in QED, it 



turns out that these effects disappear at zero momentum transfer |31[ and hence they do 
not seem to be able to give a mass to the Goldstone boson, when the latter exists. It is only 
in non-abelian gauge theories that nonperturbative effects, like those of instantons, succeed, 



through the anomalous Ward identities, in providing the Goldstone boson with a mass [p2 

Finally, a comment on the structure of the Goldstone boson is in order. This state, in 
the present mechanism of chiral symmetry breaking, is not of the same nature as that of the 
massive states. In particular, it does not result from the continuous decrease of the mass 
of a massive state down to zero, when the physical coupling constant increases, but rather 
appears abruptly as a new type of solution when the physical coupling constant exceeds a 
critical value. 
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APPENDIX: 

The wave equations of constraint theory for a fermion-antifermion system can be written 
in the form flTTl 



( 7l .pi - m x ) * = (-72-P2 + m 2 )V^, (Ala) 
(-72-P2 - m 2 ) * = (71 .pi + mi) V" * , (Alb) 

where \I/ is a sixteen-component spinor wave function of rank two and is represented as a 
4x4 matrix; the Dirac matrices 72 act on ^ from the right. The compatibility condition 
of the two equations (Al) allow one to eliminate the relative time variable and to define an 
internal three-dimensional wave function. 
After using the parametrization 

V = tanhV (A2) 

and making the change of function 

§ = (coshy) * , (A3) 

the norm of the internal three-dimensional wave function, denoted by if}, becomes (in the 
cm. frame) : 

J d 3 x Tr {^[1 + 4 7 io72oPo|^} = 2P o > ( A4 ) 

where P M is the total four-momentum of the system. 

In this representation the Todorov potential [13] takes the form (in the Feynman gauge): 



V = 7 a -72^1n(l + ^) , W = VP>. (A5) 
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Equations (Al) can be solved by first decomposing ip (the internal part of \P) on the 
basis of the matrices 1, 70, 75 and 7075 : 

= ipi + 70^2 + 75^3 + 7075^4 , (A6) 

with ipi (i = 1, ... ,4) considered as 2 x 2 matrices in the spin subspace. The relationships 
of these components with the wave function ip used throughout the text are (for the equal 
mass case and the quantum numbers s = 0, £ = 0, j = 0) : 

= ^-(Sx - S 2 ). P ^ , ^ 2 = , 

W r 

(, 2a \ ¥ .2m ( 2a \ 1/2 <p .. . 

where Si and s 2 are the spin operators of particles 1 and 2 in the 2x2 component subspace 
of the ipi's. 

In perturbation theory, a relationship can be established by means of an iterative series, 
between Eqs. (Al) and the Bethe— Salpeter equation [jr2| . In general, the potential V 
is a three-dimensional nonlocal operator in x-space, but becomes a local function when 
appropriate approximations are used. In particular, when the nonlocal operator (m 2 +p 2 ) -1 ^ 2 
is replaced by a mean value like (m 2 + < p 2 >)~ 1 ^ 2 , V becomes local in x (in the cm. frame) 
and dependent on (m 2 + < p 2 >) -1 ^ 2 . 

The Todorov potential (A5) results, however, from a slightly different approximation : 
it is a function of x and W/2, rather than of x and (m 2 + < p 2 >) -1 ^ 2 . One should then 
replace the latter quantity by W/2. It turns out that this approximation provides even 
better results, since the Todorov potential reproduces the correct spectrum to order a 4 for 
positronium and muonium [p!0| , p[| . In this approximation, at zero relative time (x° = 0), 
the relationship between the Bethe— Salpeter wave function (f>Bs and ip takes the form : 

<J) BS (x° = 0,x) = e ^™ V ^(x) , (A8) 

where V is given in Eq. (A5). The normalization constant in the right-hand side of Eq. 
(A4) takes account of this relationship. 
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FIGURES 

FIG. 1. £ = spectrum generated by Case's method for solving the KG equation (2.1) in the 
strong coupling case (aZ > 1/2), with Case's constant B = 1.56. The variable A is defined by Eq. 
(2.7). 

FIG. 2. Same as Fig. 1, with B = 1.58. The eigenvalues for all states, except the lowest one, 
are continuously continuing the eigenvalues for aZ < 1/2. 

FIG. 3. Wave function of the lowest state generated by Case's method for solving the KG 
equation (2.1), with aZ = 1 and Case's constant B = 1.58, corresponding to the dot-dashed curve 
in Fig. 2. 

FIG. 4. £ = spectrum generated by Case's method for solving the RQCD equation (3.1) in 
the strong coupling case (a > 1/2) and for Case's constant B = 1.57. The A parameter is defined 
by Eq. (3.10). All curves above and including the heavy bold curve are connected continuously 
to the spectrum obtained for a tending to 1/2 from below. As indicated in the text, all the other 
curves are merging to the origin of the axes as A — > (only the explicitly calculated parts are 
shown) . 

FIG. 5. Similar results as in Fig. 4, but for B = 1.58. The same conventions have been used. 
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